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ABSTRACT 


Time domain measurements, taken between the NAVSTAR 
6 Spacecraft Vehicle (SV) and the Vandenberg GPS 
Monitor Site (MS), by a pseudo randan noise (PIM) 
receiver, have been collected over an extended 
period of time and analyzed to estimate the long 
term frequency stability of the MVSTAR 6 onboard 
frequency standard, referenced to the Vandenberg MS 
frequency standard. 

The technique employed separates the clock offset 
from the composite signal by first applying 
corrections for equipment delays, ionospheric 
delay, tropospheric delay, earth rotation and the 
relativistic effect. The data are edited and 
smoothed using the predicted SV ephemeris to 
calculate the geometric delay. Then all available 
passes from each of the four GPS monitor stations, 
are collected at 1 -week intervals and used to 
calculate the NAVSTAR orbital elements. The 
procedure is then completed by subtracting the 
corrections and the gecmetric delay, using the 
final orbital elements, fran the conposite signal, 
thus leaving the clock offset and randan error. 

Frequency stability estimates of clock perfomance 
are then made using the clock offsets to calculate 
the Allan Variance, <^ „( ^ ), for the spacecraft 
oscillator, using sample times that vary from 1 to 
10 days. The results indicate a combined 
clock/ ephemeris frequency stability of 1.3x10 
or less, for sample times varying from one day up 
to ten days. Future work will include analysis of 
cesium standards on SV 5 and 6 as well as a 
rubidium standard on SV 5* 
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INTRODUCTION 


As part of the Navy support to the NAVSTAR GPS Clock Development 
Program, the Naval Research Laboratory (NRL) has continued (1 ) 
research and development of precise time and frequency standards. 
This paper describes the on-orbit performance evaluation of the 
NAVSTAR 6 cesium clock. 

The cesium clock in NAVSTAR 6 is the fifth one to be orbited in 
the GPS clock development program. This cesium clock was built 
under contract to the Navy by Frequency and Time Systems (FTS). 
The NAVSTAR 6 cesiua clock was activated on 26 April 1^0 (Day 
number 117, 1^0) and has been in continuous operation for more 
than one year. 

The FTS clock is a preproduction model (PM) which is designated 
as PM-11. The preproduction series of cesium beam frequency 
standards is an evolutionary development fbom the prototype model 
orbited in the Navigation Technology Satellite II (OTS-2). The 
preproduction cesium frequency standards built by FTS are 
scheduled to be placed in the NAVSTAR 5, 6, 7, and 8 spacecraft. 

The NAVSTAR 6 spacecraft is also known as SV 9* The SV 
identification is given as part of the navigation message. 

GPS SYSTM DESCRIPTION 

The NAVSTAR Global Positioning System is a Department of Defense 
(DOD) sjace-based ^stan employing a constellation of satellites 
vhich broadcast signals that are ^mchronized in both time and 
frequency. Information is given in the ephaneris message which 
can be ccmbined with measurements from four GPS satellites and 
used to calculate the user's instantaneovis position and velocity 
in all three coordinates, as well as precise time and frequency. 
Signals from the GPS satellites can be rapidly acquired and 
processed independently of all other ^stans. The precise time 
aM frequency information cai be used to provide a common, 
worldwide time grid for referencing scientific measuranents by 
laboratories all over the world. 

The NAVSTAR GPS system is comprised of four segnents: 

(1) Control Segnent 

(2) Space Segment 
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(3) User Segment 

(4) Engineering Segment 

The Space, Control, and User segments of GPS will he discussed 
with emiAiasis m factors that are iimportant to this paper. 

The GPS Control Segment consists of a Master Control Station 
(BK!S), located at Vandenherg, CA, and four reimote Monitor Sites 
(J©), located at ¥andenherg, Hawaii, Alaska, and Guam. IHiese 
fbur stations track the GPS spacecraft vehicles (SV). Data from 
these sites are collected at the MCS and processed to determine 
SV health, orbits, ani clock offsets. 

The GPS Space Segment currently consists of six NAVSTAR SVs, with 
a total of 18 satellites scheduled to be operational in the 
1986-87 time frame. These satellites are placed in (nominal) 12 
hour near-circular orbits, with occasional orbit adjust maneuvers 
Tidiich maintain a repeating gromd track for each SV. The 
configurations of the MVSTAR have been under active study since 
the original recommended constellation (2). Recent studies (3) 
have produced configurations which result in a amall tmprovament 
in the GPS coverage. 

Each GPS spacecraft broadcasts spread spectrum modulated signals 
which are precisely related to the on-board clock. The 

spacecraft navigation message (4) is also modulated onto the 

signal at precis epochs which aid in defining GPS time. 

The "precise” , or P-code modulation is generated at two 
frequenci^ in the L-band; these are designated as the L. and 
signals. The P-code signals are modulated at a rate of 10.23 
Mbps (million bits per second). The P-code modulation provides 
the cax»bility for high precision time difference measuraments, 
and is resistant to electronic countermeasures and multipath 
interferences. The P-code omploys a very long code that is reset 
once per week. The second code which is designated as the 

coarse/acquisition code, or (C/A) code, is modulated at 1.023 
Mbps and repeats every millisecond. It provides a coarse signal 
that is a factor of ten less precise than the P-code. The C/A 
signal is biphase modiilated with the P-code and may be rapidly 
acquired by all users. Each GPS SV has an atomic frequency 
standard which controls the broadcast frequency of each satellite 
to the same ncaminal value. The use of the spread spectrum 
modulation and separate codes for each GPS SV permits multiple 
access to any of the satellites that are above the user's 

horizon. 
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A GPS user would be required to have an appropriate antenna, 
receiver, processor, aM output device to receive the precise 
time and time interval signals and perform a navigation solution. 
A fully operational GPS user would select four NAVSTARs from the 
6 to 9 satellites that would he availahle in such a fashion as to 
minimize the Geometrical Dilution Of Precision (GDOP), a quantity 
[ 5 , 2 ] that relates to the navigation accuracy availahle from GPS. 
The user would acquire and lock the receiver to signals hroadeast 
from four of the SVs, and then make simultaneous measurements of 
time difference (pseudo-range) and frequency difference 
(pseudo-range-rate) between each of the SV clocks and the 
receiver clock. The user would then use the four pseudo-range 
measurements to calculate clock offset, latitiade, longitude , and 
height. The four pseudo-range-rate measurements would he used to 
calculate frequency offset and velocity in all three components, 

LONG TERM FREQUENCY STABILITY 

The GPS system is capable of providing instantaneous precise 
navigation heca^use the satellite clocks are synchronized in time 
and frequency. Therefore a fundamental measure of system 
performance is given by the long term frequency stability of each 
of the SV clocks. 

The technique that has been developed for analyzing frequency 
stability performance of orbiting clocks and frequency standards 
was developed (6) at the Naval Research Laboratory (NRL) in 1975* 

This procedure has evolved into an analytical procedure depicted 
in Figure 1. Each major ccmponent of the long term frequency 
stability analyses will be described in the developnent of this 
paper. 

FREQUENCY STABILITY MODEL 

The Allan variance was adopted by the IEEE as the reccmmended 
measure of frequency stability. Reference 7 presents a 
theoretical developnent vihidh results in a relationship between 
the expected value of the standard deviation of the frequency 
flimstuations for any finite number of data samples and the 
infinite time average of the standard deviation. Equation (1 ) 
presents the Allan variance expression for M frequency samples 
with sample period T equal to the samplir^ time, t. 

p rp .1 1 

Bn (1 ) o (M-1) 2 
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The average frequency values y-j^ are calculated from jairs <|f 
clock offsets, ^t, sej^rated by femple time, T , is given by 

E(i (2) yfc = ~ j 

I. 

The clock offset is not directly ohservable from a jsendo-rai^e 
measuronent ; other variahles must he measured or estimated. 

CIDCK DIEMEl’CE MEASUEIEMENTS 

Pseudo-range (PR) and accumulated delta pseudo-range (ADR) 
measuranents are taken between the MVSTAR SV clock and the 
clock using a spread spectrum receiver. The measurements are 
taken once every 6 seconds and then aggregated and smoothed once 
per 15 minutes. Pigiu'e 2 presents a plot of a typical 
pseudo-range signature obtained frcm a single NAVSTAR pass over a 
monitor station. Each measurement is corrected for equipment 
delay, ionospheric delay, tropospheric delay, earth rotation, and 
relativistic effects. Then the data are edited and anoothed 
using the predicted SV ephesneris to calculate the gecmetric 
delay. !Hie clock offset at the mid point of the 15 minute data 
span is estimated, using both the pseudo-range and the 
pseudo-range rate measuranents iiAiich are fitted to a ciibic 
polynoaaial with epoch at time corresponding to the mid point of 
the data. 

The pseudo-range measurements are resolved to (1/64) of a P-code 
chip, which corresponis to 1.5 nsec of time, or 46 on in range. 
Noninal values for the pseudo-range noise levels are = 1 .5m 
for the L. measurements and *^PR = 2.0m for the L 2 measuranents . 
The L. a!nd measurements are combined to correct for 
ionospheric refraction, which results in an increase to 4.55a for 
the corrected pseudo-range measurement. The accumulated delta 
pseudo-range measuranent noise levels are 0.51em for L. and 0.56 
cm for Lp. Bie^ measurements are also ccMbined to correct for 
the acetoulated pseudo-range measurements . The smoothing 
procedure uses the AER measurements to aid the PR smoothing of 
each 15 minute segnent of data. This process results in a 
anoothed pseudo-range measurement noise level of 18.5c3m. 

The equation which relates the pseudo-rar^e measuranents to the 
clock difference between the NAVSTAR SV aM the monitor site (MS) 
is given by Equation (5 ) • 

Eq (3) PR = R + c(tj^ - tgy) + c^t^ + e 
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vheve 


PR = the measured pseudo-range 

R = the slant, or geometric range, frcm the SV (at the time 
of transmission) to the I® (at the time of reception) 

c = the speed of light 

tjyjg = the MS clock time 

tgy = the SV clock time 

t. = ionospheric, tropospheric, and relativistic delay, with 
corrections for antenna and equipment delays 

e = the measurement error 


The clock difference, denoted by At, for the k.^^ measurement 
is obtained by rearranging Eq (3) into 


Eq (4) Atj^ = (tgy - tjyjg) = R/c + At^ + 0/c - RR/e 

The particular evaluation of Eq (4) that will be used in this 
paper is obtained by designating the NAVSTAR 6 by SV 9 and by 
designating the Vandenberg Monitor Site as VMS. This particular 
case of Eq (4) is given by 


Eq (5) A tj^ = (SV9 - VMS) 

S^IOOTHED ORBIT ESTIMATION 

All of the anoothed pseudo-range measurements are collected fron 
the four SPS monitor sites for one week. The Naval Surface 
Weapons Center (NSWG) then estimates a smoothed NAVSTAR orbit 
using an orbit estimation program which extensively models the 
dynamics of the satellite motion, including solar radiation 
pressure, and orbit adjust maneuvers. The NSWC post-fit 
ephemeris calculations employ the highly redxmdant set of 
range-difference (5) values, which are calculated from the 
smooldied 15 minute pseudo-range measurements. 

The purpose of the smoothed orbit estimation is to separate the 
clock and orbital ctmponents by modeling the clock as a constant. 
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(but mknown frequency) during the one week span. The model 
incorporates the feature of inclusion of an (unknown) aging rate, 
which be used for frequency standards that exhibit aging. 
The model also is capable of segjaentir^ the clock bias solution 
to allow for frequency adjustments of the B© or SV clock. 

The clock differences used for analyzing the spacecraft clock 
incorporate the smoothed orbit and the set of 15 minute 
pseudo-range measurements to calculate the clock difference at 
the time of each measurement according to Eq (4) . The clock 
differences for each MVSTAR pass are then used to estimate the 
clock differences at the time-of-closest-approach (TCA) of the 
MVSTAR SV over the monitor site. This procedure results in 

either one or two points per day. The MVSTAR orbit and the 

monitor site location determine whether one or two points per day 
will be available, Eor MVSTAR 6, one pass per day is available 
frcm the Vandenberg Monitor Site (VMS) . 

The evaluation of the clock difference at the TCA point minimizes 
the effect of the MVSTAR orbit estimation for along- the- 
satellite-track and out-of-plane errors. However this procedure 
does not reduce the effect of radial orbit errors. Hence the 
estimate of radial orbit error will be one of the factors that 
limit the accuracy of the long term frequency stal^ilty analysis. 
The effect of the radial error on the frequency stability is 
given by 


Eq (6) a ( T ) = '^^m 

^ CT 


where 

^RR = standard deviation estimate of radial component of 
orbit error 

c = the speed of li^t 
T = sample time 


MVSTAR 6 RESULTS 

The clock differences between the MVSTAR 6 cesiiM clock and the 
Vandenberg Monitor Site (1/®) clock are presented in Figure (3)* 
Each "X" ^nmbol corresponds to a single measurement obtained from 
the smoothed 15 minute pseudo-range measuranent . A total of 23 
points are plotted for this MVSTAR 6 pass over the "VMS. These 
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23 points are analyzed and a subset of these data are used to 
estimate the clock offset at the TCA point. 

The clock difference, which is denoted hy (SV9 - VMS), 
correspoMs to starting a time interval measuranent with the 
NAVSTAR 6 clock, and stopping with the VMS clock (with 
corrections for the orbit and other delays). The clock offset 
changed by approximately 30ns during the SV pass. The clock 
offset presented here represents a clock difference which ms^ be 
further processed to produce "GPS" time. The slope. ^f these 
measurements indicates a. frequency offset of -1.CBx10 between 
the NAVSTAR 6 clock and the VMS clock. Bie magnitude of this 
offset is normally what would be expected (8) after the 
correction for the relativistic clock effect. 

The clock differences for one week are presented by Figure 4- In 
this figure, each "X" symbol denotes one clock difference 
obtained from the smoothed 15 minute j®eudo-rar^e measru*anents . 
There are seven groups of "X" symbols, each one corresponding to 
a single NAVSTAR 6 pass that was observed by the Vandenberg 
Monitor Site. The slope the clock differences for this 
one-week segnent is -1.56x10 

The clock differences for the entire 100 day data span are 
presented in Figure 5. Each vertical mark corresponds to the 
clock difference evaluated at the TO A point of a NAVSTAR pass 
over the Vandenberg Monitor Site. 

Reference to Figure 5 indicates a total change in clock 
difference of about 10 Us in 100 days, or approximately 
O.lMs/day (-1 . 1 6x1 0” '^) . It also. ^di cates that a small 
frequeiKy change, on the order of 5x1 0~^, occurred between 
daj^ 120 and 150. This anall frequency shift will be further 
analyzed before computing the frequency stability. 

The frequency differences for sample times of one, three, and ten 
days are presented in Figures 6, 7, and 8. Analysis of the 
results indicates that the total frequency change between the two 
cesiWp standards during the 100 day span was on the order of 
1x10“ . Reference to Figure_.7. indicates three anall frequency 

changes on the order of 3x10 ^ occurred, with the majority of 
the frequency differences on the order of a few parts in 10(14). 
Analysis of other GPS data (not included in this ja.per) indicates 
that the Vandenberg f© cesium clock was responsible for the 
largest frequency charges analyzed. 

The frequency measurements were then used to calculate the Allan 
variance for sample times varying frcm one to ten days. An 
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interleaving (also called overlapping) data processing technique 
was used in order to obtain maximal use of the data. Ibr 
instance, with a sample time of one day and the set of clock 
differences[At. , Atp, At,, At.} two variance were calculated. 
The first varimce iSed the subset £At. ,At 2 , A t,i and the 
second variance was calculated i;^lng the subset {Atd, At_, At.} . 
Thus, the two b ( t) values have the subset { A tp, eW} ^in 
common. ^ ^ 

Frequency stability estimates, for the combined clock and 
ephemeris, are presented in Figure 9* The results are summarized 
by Table 1. 


Sample Time 
(days) 


Table 1 


Frequency Stability 
CJ ( t) Parts in 10(13) 

tX . 


1 1.1 

2 1.G 


3 

4 

5 

6 


.9 
• 9 
.9 

1.0 


7 

8 

9 

10 


1.1 

1.2 

1.3 

1.3 


The C7 ( t) may be characterized by four segnents, using a model 
given ^y 


Bq (7) ay(T) = aT^ 
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The coefficients and sample times for the combined clock and 
ephemeris frequeiKjy stability are sunmarized in Table 2. 


Table 2 


Sample Time 
(days) 

Coefficient "a” 
PP 10(15) 

Exponent ’ 

1 < T < 3 

1.10 

-0.18 

5 < < 5 

,90 

0.00 

5 < '^ < 9 

• 53 

.65 

9 < < 10 

1.50 

0.00 


Comparison of the exponents with the type of noise process 
identifiable in reference (7) in atomic frequency standards indicate 
that the segnent for sample times from three to five days is 
classified as flicker noise frequency modulation (M). The segnent 
for sample times of nine to ten days may also be classified as 
flicker noise M. The other two segments can not readily be 

classified, however the first segment has an exponent of -0.18, 
>ftiich is close to flicker noise M, and the third segment has a 
slope of 0.655, which is close to that expected frctn a randan walk 
in frequency (with an exponent of 1.0), 

The error sources that are believed to be most significant in 
limiting this analysis are 

o the use of a single monitor station frequency standard 

o the radial component of orbit error 

The first factor can be reduced by incorporating into the analysis 
multiple frequency standards at a single monitor station, or by 
analyzing the spacecraft clock performance using multiple monitor 
sites. For example, the use of the other three GPS monitor sites 
will permit the identification of frequency changes such as 
evidenced in Figure 7. A time scale could then be formed in a 
manner similar to that described in reference (8), The effect of the 
radial orbit error can be estimated by using equation (6). For the 
best fit ephemerides used in this anal^ijS, the average radial error 
was 2.1m, which corresponds to 1.4x10 for a sample time of one 
day. Additional orbit smoothing could produce better estimates for 
the orbit; however, the exponent of -1 (fran Equation (6)) is 
unchanged by additional smoothing. Ultimately, the 18cm noise level 
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obtained fran the mcwthed fifteen minute measurements yields a 
lower limit of 1.2x10 for a one day sample time. 

CONCnJSIONS 

o The frequency stabiliiy performance of the 
MVSTAR 6 cesium clock, to date, is acceptable 
and within specifications, 
o The combined (spacecraft clock, single monitor 
station clock, and ephemeris) freguency stability 
is equal to, or less than, 1.3x10 ^ for sample 
times of one to ten days. 
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Figure 1 — GPS Frequency stability analysis procedure for analyzing on-orbit clock performance 
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Figure 2— Typical NAVSTAR pseudo^range signature, 
plotted every 15 minutes 
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Figure 4— One week of clock differences between 
N AVSTAR-6 ISV9) and the GPS 
Vandenberg Morntor Site {VMS) 
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Figure 3 — Single Pass Clock difference between 
NAVSTAR-6 (SV9) and the GPS 
Vandenberg Monitor Site (VMS) 



Figure 5— 100 days of clock differences between 
NAVSTAR-6 (SV9) and the GPS 
Vandenberg Monitor Site (VMS) 
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Figure 6—100 days of frequency differences between 
NAVSTAfl-6 {SV9) and the GPS 
Vandenbeig Monitor Site <VMS) for a one 
day sample time 
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DAYS 1^1 

figure 7—100 days of frequency differences between 
NAVSTAR-6 (SV9) and the GPS 
Vandenberg Monitor Site (VMSi for a 3 
day sample time 
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Figure 8 — 100 days of frequency differences between 
N AVSTAR-6 (SV9) and the GPS 
Vandenberg Monitor Site I VMS) for a 
10 day sample time 
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Figure 9— NAVSTAR-6 Cesium Clock (PPM-11) on-orbit clock/ephemeris frequency 
stability for 1 to 10 day sample times 
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